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Abstract 

The spectral function for finite nuclei is computed within the framework 

of the Local Density Approximation, starting from nuclear matter spectral 

functions obtained with a realistic nucleon-nucleon interaction. The spectral 

function is decomposed into a single-particle part and a "correlated" part; 

the latter is treated in the local density approximation. As an application 

momentum distributions, quasi-particle strengths and overlap functions for 

valence hole states, and the light-cone momentum distribution in finite nuclei 

are computed. 
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I. INTRODUCTION 



All information on the single-particle structure of nuclei is contained in the spectral 
function S(k, E), where k and E are the momentum and the removal energy of the nucleon, 
respectively. This quantity plays a central role in a variety of high-energy nuclear reactions, 
such as in (inclusive and exclusive) quasi-free electron scattering and deep-inelastic lepton 
scattering on nuclei. With presently available many-body techniques (method of correlated 
basis functions) it can be computed quite accurately for nuclear matter starting from a 
realistic nucleon-nucleon interaction PJ . On the other hand, for finite nuclei with A > 4 it is 
much more difficult to calculate the spectral function without making severe approximations. 

In this paper we compute the spectral function for finite nuclei in the local density 
approximation (LDA) . To this end the spectral function is decomposed in terms of a single- 
particle contribution and a correlation part. The former part which varies strongly with 
mass number, is represented by the generalized mean field approximation; the latter which 
is rather insensitive to the finite size of nuclei is treated in the LDA. As an input the results 
for the nuclear matter spectral function, computed as a function of the density by Fabrocini 
et al. 0|§, are used. In order to get an idea about the sensitivity to the details of the high- 
momentum components we also use as an input a simple parametrization of the momentum 
distribution in nuclear matter given by Baldo et al. H. 

Several applications are considered. First, by integrating over the removal energy, mo- 
mentum distributions for finite nuclei with closed shells are computed. Next we consider 
three different aspects of the spectral function: (1) the quasi-particle strength for valence 
hole states and the radial shape of the overlap functions which have been measured in (e, e'p) 
reactions in several closed shell nuclei; (2) the distribution of s.p. strength at large removal 
energies where there is an important enhancement of the high-momentum components in the 
spectral function; (3) finally the light-cone momentum distribution, which plays a central 
role in the convolution approach to the EMC effect in deep-inelastic scattering on nuclei is 
computed. 
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This paper is organized as follows. In section 2 the LDA to the spectral function is 
described; section 3 deals with various applications; section 4 contains a summary and 
discussion. 



II. THE SPECTRAL FUNCTION IN THE LOCAL DENSITY APPROXIMATION 

Usually the spectral function for a finite system is expressed in terms of a (truncated) 
single-particle basis S a p(E) =< 0^5 (E — H + E )cp\0 >. Since we are interested in the 
effects induced by (short-range) nucleon-nucleon correlations on the spectral function, we 
wish to be independent of the choice and size of the s.p. basis and use the coordinate 
representation. 

For a finite nucleus the hole spectral function can be expressed as 



We use the convention €-f(A) = (-Eo(A+i) — •£'o(A-i))/2 for the Fermi energy in a finite system. 

In order to make the connection to nuclear matter it is useful to introduce the Fourier 
transform with respect to the relative distance r: 



S A (R, r, E) =<0(A)|c+(fi)<y(E -H + E 0(A) )c(r 2 )\0(A) > 



(1) 



with R = \{ri + f 2 ) and r = r\ — r<i- 



The spectral function ([!]) is normalized to the number of particles 




(2) 




(3) 



The r.h.s. of eq.(|3|) can be interpreted as the spectral function for a nucleon at position R 
with momentum k |2|||. In particular the momentum distribution in the ground state of a 



nucleus with mass A is given by 




(4) 
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For our purpose of making a LDA, it is necessary to decompose the spectral function 
and momentum distribution into a (generalized) single-particle part and a correlation part. 
The reason for this can be seen if one considers the global application of the LDA to the 
(special) case of the momentum distribution 

ni LD \k) = 1 ^JdRv(p(R),k), (5) 

where u(p, k) is the momentum density for nuclear matter at density p and with Fermi 
momentum &f(p) = (|7r 2 p) 1//3 . As noted in ||, eq.@ is divergent in the limit k — > 0. This 
problem can be circumvented by decomposing the momentum distribution into a mean-field 
part and a correlation part generated by dynamical excitations. While the former part can 
be calculated using modern mean-field theories such as density dependent Hartree-Fock, the 
latter part can then be approximated by applying the LDA to the correlation part of u(p, k) 
only, i.e. 

n A {k) = nf{k) + ^3 / dRSu^piR), k) + ^ J dR5p > (p(R), k), (6) 

with 

nT(k) = ^33 / dRdfe-^nTiR^, 
5u<(p, k) = (z/(p, k) - l)6(k F - k) 

6v > (p,k) = u{p,k)9(k-k F ), (7) 
and the mean-field density given by the sum over the occupied orbitals 

nT(R,r) = X>r(^) = E«n)<Mr 2 ). (8) 

h=l h=l 

A similar divergence problem, which occurs if the LDA is applied to the full spectral 
function, can be treated in an analogous way. In this case there are several ways to decompose 
S in different parts. Let us first consider only the mean-field part of the interaction felt by 
a nucleon in the nucleus, in which case the spectral function has the simple form: 

S^(R, f, E) = J2 nf F (R, f)8(E + e h ), (9) 

h=l 
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In addition to the hole orbitals appearing in eqs.([8|-|j|), the mean-field also defines the nor- 
mally empty particle orbitals <p p . As for deviations from the mean-field approximation it is 
convenient to distinguish the following two dynamical effects of the residual interaction: 

1. Depletion and fragmentation of the hole strength. 

Nucleon-nucleon correlations lead to both depletion and fragmentation of the single 
hole strength in eq.@. We assume that this depletion and fragmentation of the 
strength related to the hole orbital h only depend on the single hole energy, e^; they are 
taken equal to the depletion and fragmentation of the spectral function S nm (p, k, E) in 
nuclear matter at the local density. Consequently the 5-function in the single-particle 
part of the correlated spectral function in eq. @ is replaced by the distribution function 

8{E + e h ) -> S nm (p(R), k(e h ), E). (10) 

This ensures a realistic behaviour of the strength distribution at large removal energies; 
it is well known that Lorentzian distributions become inadequate away from the quasi- 
particle peak and cannot be used for the evaluation of e.g. mean removal energies. The 
momentum k(eh) < hp corresponding to is chosen in such a way that the s.p. energy 
(relative to the Fermi energy) ef(A) — e h corresponds to the position of the quasi-hole 
peak in the nuclear matter spectral function. The fact that the finite volume of the 
nucleus imposes discrete energies and bound single-particle wave functions 4>h (of well- 
defined angular momentum) is taken into account by keeping the mean-field density 
nJf F (-R, r) in eq.@. The single-particle part is dominated by the mean field, and the 
^-dependence is mainly determined by the hole orbitals. 

2. Partial occupancy of the normally empty states. 

The scattering of nucleon-nucleon pairs into states above the Fermi momentum leads 
to both a depletion of the occupied Fermi sea and to the presence of high-momentum 
components in the spectral function at large removal energies. Since this effect is 
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expected to depend only weakly on the shape of the orbitals P in the finite nucleus 
the correlation part of Sa is taken equal to the spectral function in nuclear matter for 
k > kp at the local density p(R): 

ST(R, k, E) = ^S nm (p(R), k, E)6(k - k F (R)). (11) 

We emphasize again that this prescription cannot be extended to the whole spectral 
function (k < kp) since an integration over R would lead to a singular result for the 
momentum distribution at k — (see eq.(^)). The separation of the spectral function 
into a single-particle part and a correlation part is thus quite natural, as the low- 
momentum part of the spectral function depends crucially on the finiteness of the 
system. 

In summary, we use the following decompositon of the hole spectral function 
S A (R, k, E) = S S A P (R, k, E) + S c A or (R, k, E) 

with 

S S A P (R, k,E) = J2 nf Y {R, k)S am {p{R), k{e h ), E) 

h=l 

and S c A or (R, k, E) given by eq.(0). 

III. APPLICATIONS 

A. Nuclear matter input 

For the nuclear matter spectral functions we used the results obtained by Fabrocini et 
al. [0,|J for densities p/p = 0.25, 0.5, 0.75 and 1.0, with p = 0.16 fm -3 the nuclear matter 
density at equilibrium. The spectral functions were calculated with correlated basis function 
(CBF) theory using the Urbana Vu NN interaction with addition of a three- nucleon force 
0. The calculations involve a variational determination of the correlated groundstate and 
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(12) 



(13) 



of the lh and 2hlp states. Furthermore perturbative corrections (2p2h) to the groundstate 
and (2hlp) to the lh states are added. 

As noted in ref. |3[] a peculiar feature of the nuclear matter results is that the momentum 
density does not approach the non-interacting one in the limit of zero density. For example, 
one finds a rather constant depletion, e.g. u(p, k « 0) = 0.86, 0.87, 0.86 for p = p , |p , and 
|po; respectively, and similar results for the discontinuity Z(p) of u(p, k) at the Fermi surface. 
This behavior can be attributed to the fact that in the variational approach the attraction 
of the NN interaction leads to bound pairs and thus to clustering for low densities. One may 
argue that this effect is an artefact of the variational approach and is not a realistic density 
dependence to be used in the LDA for finite nuclei. Therefore we assumed that 8v(p, k) — > 
for p — * by a smooth extrapolation to zero for densities lower than jp . One should note 
that the final results are quite insensitive to this procedure because the weight of small p 
values is quite small. 

To investigate the sensitivity to the nuclear matter input, it is of interest to compare the 
results with those for simple parametrizations of the nuclear matter momentum distribution 
as a function of the density. Here we use the parametrization recently given by Baldo et al. 



Su(p,k) 



-0.21 - 0.13k - 0.19(1 - k) ln(l - k), k < 1 

0.21 + 0.3 arctan x + 0.82x \nx, 1 < k < 2, k< 2fm _1 (14) 



% exp (-1.6&), 2 < k < 4.5ml- 1 



7.0 

where x — (« — 1)/ (k 2 + 1) and k = k/kp. This is an approximate parametrization valid for 
k-p > lfm , in the sense that for k < 2fm _1 the result depends on k only. Since we require 
that 5v(p, k) — > for lower densities, we multiplied the above expression for Su(p, k) with k 2 F 
for kp < lfm -1 . We also properly renormalized the momentum distribution for all densities 
considered. The parametrization of eq. (0) provided a good fit to nuclear matter momentum 
distributions with 1 < kp < 1.75 fm _1 calculated to second order in the Brueckner reaction 
matrix and using a separable representation of the Paris interaction. 
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In Figure 1 we compare u(k) from the two prescriptions for nuclear matter at two den- 
sities. Note that the former has a higher depletion and somewhat stronger correlations. 

This is also reflected in Table 1, where we list the kinetic energies obtained with the 
spectral functions of |2j and with the parametrization in Q], as a function of the Fermi 
momentum in nuclear matter. For the kinetic energies obtained with the spectral functions 
of |2| we assumed an exponential extrapolation of the momentum distribution for large 
momenta. We note that the total kinetic energy in the interacting system is roughly twice 
the free value given by T free = y^. 

B. Momentum distributions in finite nuclei 

The momentum distributions for finite nuclei are computed by integrating eq.(|12|) over 
R and E 

n A {k) = nf(k) + J dR^nfViR^SvMR),^)) + ^ j dR6 V> (p(R),k), (15) 

for the spherical closed-shell nuclei with A =16, 40, 48, 90 and 208. We checked that the 
normalization @ was in all cases fulfilled to better than 1%. The mean-field densities were 
taken from @. For comparison we also used the more conventional form of the LDA of 
eq.@ for the momentum density. The two expressions (eq.@ and eq.(|i~5"D) differ in the 
second term, which represents the depletion of the mean-field momentum distribution at 
small momenta. In Figure 2 the resulting momentum densities are compared for 16 and 
208 Pb. As expected, the high-momentum components (k > 1.6 fin -1 ) are identical for both 
cases. In the intermediate range (lfm -1 < k < 1.6 fm" 1 ) the conventional LDA gives rise to 
a kink in the momentum distribution. The kink is due to the second term in eq.(^), which 
(with fcf max the Fermi momentum corresponding to the largest density in the nucleus) has 
a very steep behaviour for k f> &F m ax an d is strictly zero for k > kp max . This unphysical 
feature is removed in the present approach, since the second term in eq.(|l~5"D has a smooth 
behaviour for all k. The present treatment also predicts for 16 somewhat larger riA{k) for 



S 



k ~ than the conventional LDA, because eq.(|15"D takes into account that the reduction of 
the mean-field distribution for k ~ comes only from the depletion of the lsl/2 (I = 0) 
orbitals, which are deeply bound and therefore less depleted than the lp orbitals. The 
conventional LDA, on the other hand, can only take into account an average depletion. For 
16 we also included in Figure 2 the result of a Variational Monte-Carlo calculation by 
Pieper et al. ||. As was noted in H, there is good agreement up to k = 2fm _1 between the 
LDA and the VMC approach. For k > 2fm _1 the difference can probably be ascribed to 
the fact that the VMC calculation was carried out using the Argonne NN interaction, which 
has a stronger tensor force than the Urbana interaction ||. 

In Figure 3 the LDA result for the A-dependence of the momentum density is shown. 
Although with increasing A there is a clear tendency towards the nuclear matter momentum 
density, finite-size effects are seen to remain important (even for 208 Pb), especially at small 
k. 

The resulting A-dependence of the kinetic energy T, AT = T — T MF , and the mean 
removal energy < E > is given in Table 2. It is seen that both the MF and the correlation 
part of the kinetic energy amount to approximately 17 MeV, and increase with A. The 
difference between AT in the present approach and in the conventional LDA is small (at 
most 3%), whereas the parametrization of the nuclear matter momentum density by Baldo 
et al. leads also in finite nuclei to sizeably larger AT. The mean removal energies < E > 
in Table 2 were calculated using the mean-field (|) and correlated spectral function (|12"|), 



and have been corrected for the mean-field rearrangement energy. The inclusion of nucleon- 
nucleon correlations increases < E >, from 33 MeV in the mean-field approximation to about 
55 MeV. Note that the effects of the Coulomb interaction and the proton-neutron asymmetry 
are only taken into account via the mean-field single-particle energies. For completeness we 
have also compared in Table 2 the experimental binding energies per particle with 
the binding energies following from the Koltun sum rule 

Ea/A=U^T-<E>), (16) 



2 \A 



which is valid if only two-body forces are present. We find about 2.4 MeV too much binding. 
This discrepancy is not too worrisome since we neglected the effect of three-body forces, and 
(part of) the Coulomb interaction and proton-neutron asymmetry effects. 



C. Spectroscopic factor and radial shape of the quasi- hole wavefunctions 

1. The LDA approach 

Of special interest is the spectral function for the least bound orbitals just below the 
Fermi energy. For values of the removal energy near the Fermi energy e p the nuclear matter 
spectral function vanishes except for momenta slightly below kp. As a consequence the 
correlation part S A OT in eq.flTTD vanishes and the spectral function in LDA can be expressed 
as 

S A (R, E « e F{A) ) = <P* h (R)MR)Z(p(R))S(E + e h ). (17) 

h 



By comparing eq. ( |T7j ) with the general form of the spectral function for the top shells in the 
(A — 1) nucleus 

S A (R, E) = J2 I < <~V{R)\^a > \ 2 5(E + e h )=J2 \MR)\ 2 S(E + e h ) 
h h 

one sees that in the LDA the overlap functions iph{R) for the low- lying quasi-particle states 

can be expressed as 



tPh DA {R) = sJZ{p{R))<j> h {R). (18) 



This result (|18|) has the correct asymptotic behaviour of the overlap function which is deter- 



mined by the separation energy of each orbital |J. It has also been noted that for the case 
of droplets of 3 He atoms eq.(^) indeed provides a good approximation to the exact overlap 
functions | fTU| . 

Recent (e, e'p) and (7,p) experiments (after correction for MEC) have probed the spectral 
function for the top shells at large missing momenta. Therefore it is of interest to examine 
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to what extent the overlap functions for valence hole states are affected by high-momentum 
components due to short-range NN interactions. One sees from eq. flIB|) that in the LDA the 
only modification of the mean-field wavefunction in coordinate space can be expressed in 
terms of a modulation with a density dependent factor that reduces the s.p. strength. 

As a typical example Figure 4 shows that the net effect on the radial shape of the lpi/2 
overlap function in 16 O is a slight depletion of the interior and an enhancement of the surface 
and tail region. In momentum space this leads to only small global deviations from a typical 
mean-field wavefunction. This is in agreement with results obtained by Miither and Dickhoff 
||12|| , who computed the overlap function for the lpi/2 orbital in 16 using Green function 
perturbation theory with a G-matrix interaction derived from a realistic NN interaction. 

The spectroscopic strength of a quasi-particle state h is in this approach given by: 

Z LDA = J dRZ(p(R))\MR)\ 2 - (19) 

From Figure 5 one sees that in the LDA the spectroscopic factors (averaged over the valence 
hole states) decrease slowly with increasing A. This is related to the increase of the mean 



density with A. In comparison the data of ref. JXT[] also show little A-dependence, but the 



experimental spectroscopic factors are about 0.1 smaller than the LDA prediction. 



2. Inclusion of surface effects 

One effect not taken into account in a LDA treatment of correlations is the occurrence 
of surface degrees of freedom in finite nuclei. In particular collective low-lying excitations 
in the core nucleus (surface vibrations) are known to couple strongly to the valence s.p. 
states, and may affect the shape of the radial overlap function. Namely the valence hole 
state in the residual (A-l)-nucleus has one-phonon/one-hole components which can have a 
non-vanishing single-particle overlap with two-phonon components in the ground state of 
the target nucleus. This leads to admixtures in the overlap function of s.p. states which, in 
contrast to the hole mean-field wave functions, do have sizable components in the momentum 
range 1.5 < k < 2.5 fm -1 . 

11 



Following refs. [14 , 15] we have included the effect of surface vibrations on the overlap 
functions in 208 Pb through the energy dependence of the mean field. In general both the 
non-locality and the genuine energy dependence can be included by the introduction of a 
total effective mass m*(R) as the product of a fc-mass (rrik) and an E-mass (mg) 

m*(R) m k (R)m E (R) 



m mm 



(20) 



Since we assume that the non-locality (giving rise to m k ) is already included in the Skyrme 
HF wavefunctions and that all volume contributions to mg are already incorporated in 
the LDA, we only take into account explicitly the surface peaked .E-mass m^(i?)/m = 
(1 + fisdg/dR). Here the central part of the mean-field potential g(R) and the parameter 
Ps, which reflects the coupling strength of the surface modes to the s.p. states, were taken 



from | 15| . We find that he resulting overlap function, ip^ (R), which includes the effect of 



surface vibrations, can in good approximation be expressed by an extra modulation factor 



in eq.(|i"8|) 



^y(R) « x r^^iJz{ P {R))4> h {R). (2i) 



lm E (R) 



We have checked e.g. that for the valence proton s.p. states in 208 Pb the quasi-particle 
wavefunctions in are well approximated by the standard Woods-Saxon wavefunctions 
multiplied by the effective mass according to (0), although the surface peaking is somewhat 
overestimated. We note that this prescription was also proposed in ref. ||13|| . 

Figure 6 shows that in momentum space the momentum distributions generated by 
eq.(|2lD are enhanced compared to overlap functions without surface effects, by a factor of 
10-100 for momenta in the range 1.5 < k < 2.5 fm -1 . There is indeed recent evidence 
from an (e, e'p) experiment on 208 Pb [|n| for a systematic enhancement of the momentum 
distribution in this range, as compared to mean-field valence hole wavefunctions. However, in 
a detailed analysis of the experimental data it should be taken into account that the inclusion 
of surface effects through a surface-peaked modulation factor changes e.g. the r.m.s. radius 
of the overlap function, as well as the momentum distribution at smaller momenta. 
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The extra reduction Z sv of the quasi-particle strength that results from the coupling 
to surface vibrations can be calculated from the expectation value of (m^fi^/m)" 1 with 
respect to the overlap function ip^ v . For the total quasi-particle strength Z = Z sv .Z LDA 
we obtain an average of Z = 0.50 for the valence hole shells in 208 Pb. This result is similar 
to the one obtained in [IB|, but is 0.15 lower than the data in [IT]. This is probably due to 
double counting when adding the surface effects to the nuclear matter results. 



D. Global distribution of single-particle strength 

The energy distribution of total proton hole strength 

S A {E) = JdkS A (k,E) (22) 

is shown in Figure 7 for 208 Pb, split in s.p. and correlation parts. The s.p. part is dominant up 
to 50 MeV removal energy. It represents the quasi-hole strength corresponding to knockout 
from the various proton hole shells. Beyond the energy region of the quasi-hole excitations 
the s.p. part of the strength dies out quickly, and the correlation part of the strength, which 
extends to very high removal energies, becomes dominant. 

Also of interest is the difference in the momentum distribution between small and large 
values of the removal energy. Figure 8 clearly illustrates that high-momentum components 
are correlated with large removal energies. Qualitatively this is easily understood, since 
in order to remove a nucleon with high-momentum in the groundstate one has to break a 
correlated pair. The remaining nucleon in the pair has roughly opposite momentum, its 
kinetic energy leading to high excitation energies in the residual system. 



E. Light-cone momentum distribution 

The light-cone momentum distribution plays a central role in the description of high- 
energy reactions on nuclear targets, such as deep-inelastic scattering. For example, in the 
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simple convolution approach (use of the impulse approximation, and neglecting off-shell ef- 
fects) to the EMC effect the structure function F 2 (x) for the A-body target can be expressed 
as [0 

poo 

F 2 A (x) = / dyf A (y)F 2 N (x/y) 

J X 

where the nuclear structure information is contained in the light-cone momentum distribu- 
tion 

f A (y) = J dEd 3 kS A (k,E)y5(y-^-^), (23) 

where the factor y represents the flux factor. For values of x < 0.6 one does not need the 
full spectral function S A (k,E), and it is sufficient to expand F 2 (x/y) around y = 1 

F 2 {x/y) = F 2 {x) -{y- l)xF' 2 {x) + l -{y - l) 2 {2xF 2 {x) + x 2 F 2 (x)), (24) 

and therefore only the lowest moments of f(y) are needed. The lowest moments of f(y) (up 
to order can easily be expressed in terms of T and < E > 

<E>+\T 2 2T 

<y-l>= 3_, < 1 - y) 2 >= — . 

m im 

with the result 

F 2 \x)/A = F 2 N (x) - ^^xFf(x) + ^Ff{x). 

In the earlier analysis of the EMC effect for finite nuclei mostly the MF result has been used 
for < E > and T ||18|| , which can explain only 30% of the observed reduction at intermediate 
x. The inclusion of correlations enhances the effect in the ratio Ra(x) = F 2 (x)/AF 2 (x) by 
about a factor 2, but still cannot fully explain the observed EMC ratio as a function of x 

0- 

It is also of interest to examine whether the mass dependence predicted by the convolution 
method agrees with the observed one. To this end it is convenient to parametrize the 
A dependence of T and < E > in terms of a 'volume' and a 'surface' term, i.e. T = 
T vol — T surf A -1 / 3 , and similarly for < E >, hence Ra(x) = ci(x) + c 2 (x) A -1 / 3 . As an example, 
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using the best fit values to the LDA results for T vol (=39.38MeV) and T surf (=18.87MeV), 
and using the Koltun sum rule for the mean removal energies, one finds for x = 0.60 that 
R A (x) = 0.83 + 0.18A-V3. 

Recently a new analysis of the A-dependence of deep-inelastic electron scattering at 
SLAC was presented |]20|| . In that analysis the EMC ratio for the range of A considered 
(2 < A < 197) was parametrized as a power law: R A (x) = C(x)A a ^ with C(x) « 1.0; 
e.g. at x = 0.60 the best fit is obtained for a(x) = —0.0346 ± 0.002. We note that this 
parametrization has a rather unphysical limit for A — > oo. If we use the more physical 
parametrization above we obtain c\ = 0.82, C2 = 0.23 which gives an equally good fit to the 
data (see Figure 9) and moreover is in good agreement with our prediction. This indicates 
that the observed A dependence of the EMC effect is consistent with a purely single nucleon 
effect. On the other hand if two-body effects would be important a more general mass 
dependence would be expected. 

IV. SUMMARY AND CONCLUSIONS 

In this paper we have proposed a method to obtain the single nucleon spectral function 
for finite nuclei from that of nuclear matter by applying the local density approximation. 
Recently a similar method for the construction of the spectral function in finite nuclei was 
proposed @ (with the same nuclear matter input) and applied to inclusive electron scattering 
at large q. The decomposition of the spectral function in || is different from ours; e.g. at the 
variational level the s.p. part only contains the direct g.s. — > lh contributions, whereas the 
g.s. — > 2hlp background contributions are part of the correlation part for both k < kp and 
k > kp. On the other hand in our approach the correlation part contains the full nuclear 
matter spectral function only for k > kp, whereas the background contribution for k < kp 
is incorporated in the generalized single-particle part. However, the resulting momentum 
distributions for 16 O from |2j and from the present work are almost identical. 

A point of uncertainty in the LDA treatment of nuclear correlations is the extrapolation 
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of the nuclear matter spectral functions in the limit of zero density, but this is expected to 
be of little influence on the final results. 

We find that correlations lead to an appreciable depletion of quasi-particle strength Z 
in nuclei; in r-space the valence overlap orbital is somewhat more reduced in the interior 
than in the exterior region, but we find that the momentum distribution of the valence 



orbit overlap is hardly affected. This agrees with the conclusion in ref. |L2[]. However, 
surface effects may lead, for intermediate momenta, to a considerable enhancement of the 
momentum distribution for valence hole states. 
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TABLES 

TABLE I. Kinetic energy [MeV] per particle as a function of Fermi momentum in nuclear 
matter, calculated with the momentum density from Baldo et al. Q, and from Fabrocini et al. ||]|| 



ftp (fin -1 ) 


1.13 


1.23 


1.33 


Tf re e 


15.88 


18.82 


22.00 


Baldo 


37.62 


42.93 


48.31 


Fabrocini 


33.25 


36.99 


40.77 



TABLE II. The A-dependence of kinetic and mean removal energies [MeV] per particle. T 
and < E y MF refer to the mean-field approximation. AT a ^ were calculated with the nuclear 
matter input from Fabrocini et al. HQ, using ( a ^) the present model for the spectral function and 

the conventional LDA. For AT C ' the nuclear matter input from Baldo et al. j| was used, in the 
conventional LDA. The mean removal energy < E > was calculated with the present LDA for the 
spectral function and includes the mean-field rearrangement energy. Ea/A is the binding energy 
per particle, calculated using the Koltun sumrule with < E > and T a ' . 



A 


16 


40 


48 


90 


208 


rpMF 


15.4 


16.5 


17.7 


17.9 


18.6 




16.3 


16.9 


17.0 


17.3 


17.5 


AT b) 


15.8 


16.5 


16.7 


17.0 


17.3 


AT C ' 


18.5 


20.1 


20.7 


21.4 


22.2 


<E> MF 


30.4 


33.2 


34.8 


35.1 


34.2 


<E> 


50.4 


54.0 


55.9 


56.5 


56.6 


E A /A 


-10.6 


-11.0 


-11.1 


-11.0 


-10.5 


E A /A(exp.) 


-8.0 


-8.6 


-8.7 


-8.7 


-7.9 
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FIGURES 

FIG. 1. Nuclear matter momentum distributions (normalized to |vr/cp) at densities p/po = .5 
and 1, according to the spectral function of Fabrocini et al. (2],|j (full line) and with the parametriza- 
tion of Baldo et al. (dashed line). 

FIG. 2. Momentum density (normalized to unity) in 16 and 208 Pb. Short-dashed line: 
mean-field result. Full line: present LDA treatment. Long-dashed line: conventional LDA. The 
dots in the 16 O plot were taken from (8|. 

FIG. 3. Momentum density (normalized to unity) for finite nuclei (dashed lines) and nuclear 
matter at equilbrium density (full line), calculated in the present LDA treatment, with the nuclear 
matter input taken from f2|,f|. 

FIG. 4. Overlap function (normalized to unity) for the 16 0-lpl/2 state in r-space and g-space. 
Dashed line: mean-field approximation. Full line: including correlations in the LDA according to 
eq.(PJ). 



FIG. 5. Average spectroscopic factor for the valence hole states, as a function of mass number. 
The LDA prediction of the present work (squares) is compared with data (circles) taken from a 
compilation of (e, e'p) measurements in |pdj| . 

FIG. 6. Overlap function (normalized to unity) in (/-space for valence proton hole states in 
208 Pb. Short-dashed line: mean-field approximation. Full line: including correlations in the LDA 
according to eq. fll~8|) . Full line with dots: including additional surface effects according to eq.(^). 

FIG. 7. Proton spectral function for 208 Pb, integrated over all momenta, as a function of 
removal energy. Short-dashed line: single-particle part. Long-dashed line: correlation part. Full 
line: total. 
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FIG. 8. Proton spectral function for 208 Pb, integrated over various regions (0-50, 50-100, 
100-150 and 150-200 MeV) of removal energy, as a function of momentum. 

FIG. 9. Mass dependence of the EMC ratio Ra(x), for x = 0.6. The SLAC data and best power 
law fit (long-dashed line) were taken from |2(|. Short-dashed line: best fit with the parametrization 
in terms of a volume and surface contribution (see text). Full line: LDA prediction for this 
parametrization. 
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